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Abstract 

The gauge invariant quark Green's function with a path-ordered phase factor along a 
straight-line is studied in two-dimensional QCD in the large- N c limit by means of an exact 
integrodifferential equation. Its spectral functions are analytically determined. They are 
infra-red finite and lie on the positive real axis of the complex plane of the momentum 
squared variable, corresponding to momenta in the forward light cone. Their singularities 
are represented by an infinite number of threshold type branch points with power-law 
—3/2, starting at positive mass values, characterized by an integer number n and increas- 
ing with n. The analytic expression of the Green's function for all momenta is presented. 
The appearance of strong threshold singularities is suggestive of the fact that quarks could 
not be observed as asymptotic states. 

PACS numbers: 12.38.Aw, 12.38.Lg. 
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1 Introduction 



Conventional quark Green's functions suffer from lack of gauge invariance. In particular, 
no physical conclusion can be drawn from any of their individual properties. Their role 
is rather limited to lead, together with other Green's functions, to the evaluation of 
gauge invariant observable quantities, like bound state masses, current matrix elements, 
hadronic form factors, scattering amplitudes, etc. This is not the case of the gauge 
invariant quark Green's functions (GIQGF), constructed with the aid of the gluon field 
path-ordered phase factors [U [2] . Their properties should provide a more stable basis for 
the investigation of the physical characteristics of the underlying quark and gluon fields. 
Nevertheless, because of the mathematical complexity of the corresponding composite 
objects, not much has been known up to now about them. 

Recently, the present author could obtain, using a representation for the quark propa- 
gator in an external gluon field that generalizes the one introduced in the nonrelativistic 
case [3], exact integrodifferential equations for the the two-point GIQGFs (2PGIQGF), in 
which the paths of the phase factors are of the skew-polygonal type [I] . It turns out that 
such GIQGFs form a closed set with respect to the above equations. The kernels of these 
equations are represented by functional derivatives of Wilson loop averages (5j El [TJ El E] . 
Designating by Sm(x,x') = S(x,x') the 2PGIQGF with a path made of one straight line 
joining the quark to the antiquark and by Sua(x, x'\ y n -i, ■ ■ ■ ,Ui) (n >2) the 2PGIQGF 
with a path made of n segments and n — 1 junction points yi,y 2 , ■ ■ ■ ,Hn-i between the 
segments, the equation satisfied by 5 has the following structure: 

{i^.dix) -m)S(x,x') =id 4 (x-x') + iY {K 2ft (x', x, y x ) S( 2 )(yi, x) 

oo 

+ J2 K t^(x',x,y u . . .,yi-i) S {i) (yi_ u x';x,y u . . .,j/i_ 3 )}. (1.1) 

i=3 

[m is the quark mass parameter. Integrations on the intermediate variables y are implicit.] 
The kernels K n (n = 2, . . . , oo) contain logarithms of Wilson loop averages along the 
contour of the (n + l)-sided skew-polygon with at most n functional derivatives along 
the contour and (n — 1) GIQGFs S along the segments of the contour. The 2PGIQGFs 
S( n ) { n > 2) are themselves related to the lowest-order 2PGIQGF S through functional 
relations involving mainly the Wilson loop average along the (n + l)-sided skew-polygon. 
It is then conceivable that, by eliminating the S( n )S in favor of S, one might reach a closed 
form equation for S. 

What kind of information could provide us the analysis of the 2PGIQGF S7 The 
status of that Green's function is rather intricate. If colored objects are confined and 
if singularities of Green's functions result from insertions of complete sets of physical 
states, which are here composed of hadronic states (color singlets), then it is not possible 
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to cut the line joining the quark to the antiquark and saturate it with physical states, 
the contribution of the latter being identically zero. This situation would suggest that 
the above Green's function is free of singularities. Nevertheless, the equation that it 
satisfies and which is derived from the QCD lagrangian signals the presence of singularities, 
generated by the free quark propagator. This paradoxical situation is solved only with 
the acceptance that the building blocks of the theory, which are the quark and gluon 
fields, continue forming a complete set of states with positive energies and could be used 
for any saturation scheme of intermediate states. It is then up to the theory to indicate 
us through the solutions of the equations of motion how these singularities combine to 
form the complete contributions. For sectors of colored states, we should expect to find 
singularities, whenever they are present, which would be stronger than the observable 
ones. Therefore, the determination of the 2PGIQGF would provide us a direct test of the 
confinement mechanism in the colored sector of quarks. 

Saturating the intermediate states of the 2PGIQGF S with quarks and gluons with 
positive energies, one arrives at a generalized form of the Kallen-Lehmann representation 



[S(p) is the Fourier transform of S(x,x') = S(x — x').} The increasing powers of the 
denominators come from the increasing number of gluon fields in the series expansion 
of the phase factor in terms of the coupling constant. Formally, one might reduce, by 
integrations by parts, all the denominators to a single one, but the possible presence of 
strong threshold singularities might stop the operation at some power. In the latter case, 
another alternative would be the extraction out of the integral of derivation operators. 

There might also exist other extreme cases. For instance, the singularities of S, as 
provided by the solution of Eq. (II. ip . might lie in unphysical regions, like the complex 
plane or the negative axis of p 2 . Or, the summation of the series (II. 2p might lead to 
new singularities in unphysical regions. We shall, however, continue exploring the opti- 
mistic scenario, according to which the series of singularities (II. 2p are the only type of 
singularities of S such that their sum remains on the physical timelike axis of p 2 . 

This paper is devoted to the resolution of Eq. (II. ip in two-dimensional QCD in the 
large- N c limit [T0| [TT| IT2] . That theory displays the main features of confinement of 
the four-dimensional case with the additional simplification that asymptotic freedom is 
realized here in a rather trivial way, since the theory is superrenormalizable. Furthermore, 
many aspects of two-dimensional theories are often exactly solvable and in this case one 
might hope to have an unambiguous check of the physical implications of Eq. (II. ip for 
the spectral properties of the 2PGIQGF. 

The main results of the present investigation can be summarized as follows. The series 
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of spectral functions of the vector and scalar types of Eq. (I1.2p can actually be grouped 
into two single spectral functions. The singularities of the latter can be analytically 
determined. They are represented, by an infinite number of threshold type branch points 
with power-law —3/2, starting at positive mass values characterized by an integer number 
n and increasing with n. This situation is qualitatively unchanged in the limit of massless 
quarks. The strong thereshold singularities, together with the absence of pole terms, are 
suggestive of the fact that quarks could not be observed as physical asymptotic states. 

The plan of the paper is the following. In Sec. [21 we consider two-dimensional QCD 
in the large- iV c limit. In Sec. El the Green's function is considered in x-space and its 
properties are determined. In Sec. SI the instantaneous limit of the Green's function in 
momentum space is taken. In Sec. [51 the analyticity properties of the Green's function 
in the instantaneous limit are studied and the spectral functions are reconstructed. Sec. 
El is devoted to a brief application to the calculation of the quark condensate. Conclusion 
follows in Sec. [71 An appendix is devoted to the study of the properties of the threshold 
masses and to their determination. 



2 Two-dimensional QCD in the large-A^ limit 

On technical grounds, considerable simplification occurs in Eq. (11 -ip in two-dimensional 
QCD at large N c . Wilson loop averages are now exponential functionals of the areas 
of the surfaces lying inside the contours [131 E]- First, for simple contours, made of 
junctions of simpler ones, the Wilson loop average factorizes into the product of the 
averages over the simpler contours. This feature leads to the disappearance in the kernels 
of Eq. (II. ip of nested-type diagrams. Second, crossed type diagrams as well as quark loop 
contributions disappear. Third, the second-order functional derivative of the logarithm 
of the Wilson loop average is a (two-dimensional) delta-function. As a consequence, since 
each derivation acts on a different segment, all functional derivatives of higher-order than 
the second vanish. The series of kernels in the right-hand side of Eq. (II. ip then reduces 
to its first term. Equation (II. ip becomes: 

(i r d (x) - m)S{x - x') = i5\x - x') - i^f ( ^!^f } ) S{x - y) i a S {2) (y, x>; x), (2.1) 

where F3 is the logarithm of the Wilson loop average along the triangle xyx'. The func- 
tional derivatives acting on F 3 represent integrated derivatives on the segments xx' and 
yx', x' remaining fixed [I]; one has: 
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S 2 F 3 (x',x,y) 
5x^~5y a+ 



x f 1 dX dX' X X'5 2 (X'(y - x') + X(x' - x)), 



(2.2) 



where a is the string tension. The 5-function forces the point y to lie on the line x'x. 
In the present theory, the 2PGIQGF Sty is related in a simple way to S: S( 2 )(y,x'] x) = 
e F 3 (x',x,y)g^y _ x >y Integrating with respect to y and one of the As with an integration by 
parts, one obtains the following (exact) equation: 



(i^.d - m)S(x) = iS 2 (x) - ^{g^gup - g^gv a )x v x p 



x 



dX X 2 S{(1 - X)x)^ a S{Xx) + / d£ S{{\ - £)x)-? a S(£x) 



(2.3) 



The first integral represents contributions when y is lying inside the segment x'x, while 
the second one contributions from outside the segment x'x in the direction x'x. 



3 The quark Green's function 

Let us first observe that attempting to develop Eq. (12. 3p into an iterative series with 
respect to the free propagator would lead to a series of terms with increasing singularities 
at the quark mass m. Replacing in the interaction dependent part of Eq. (12. 3p S by 
the free propagator Sq, produces a term that behaves at the quark mass threshold as 
(p 2 — m 2 )~ 5 / 2 , which, in turn leads to stronger singularities at higher orders. Therefore, 
the existence of a nonperturbative solution, with possibly stable properties concerning its 
singularities, would represent the only issue to have access to the structure of the complete 
solution. It turns out that Eq. (12.31) does have such a solution satisfying the required 
boundary conditions and expressible in analytic form. A complete understanding of the 
spectral properties of the Green's function S requires an analysis in momentume space. 
We defer that study to Sees. H] and [5j In this section, we shall display the solution and 
shall check in x-space that it satisfies Eq. (12. 3p . 

Decomposing the momentum space 2PGIQGF S(p) along vector and scalar invariants 
Fi and F , 

S(p)= 1 .pF 1 (p 2 ) + F (p 2 ) } (3.1) 
the solutions for Fi and F , for complex p 2 , are: 

OO 

^) = -^K (Ml _^ (32) 

(3-3) 
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where the coefficients b n are parametrized as 



(3.4) 



M n + m n + (-l) n m 

The mass parameters M n have real positive values, greater than the free quark mass m, 
increasing with n, while the parameters m n generally have real positive values decreasing, 
for a fixed parity of the index n, with n. They are defined through an infinite set of 
algebraic equations which will be displayed below. Their asymptotic values, for large n, 
such that n ^> m 2 / (cnr), are: 

Ml ~ airn, m n ~ --jrp. (3.5) 

7T M n 

The functions (M 2 — p 2 ) 3 / 2 are defined with cuts starting from their branch points 
and going to +00 on the real axis; they are real below their branch points down to —00 
on the real axis. 

From Eqs. (13 .2p and (I3.3p we deduce that the singularities of the functions Fx and 
Fq are located on the positive real axis of p 2 and represented by an infinite number of 
branch points or thresholds with stronger singularities than simple poles. That feature is 
an indication that the resulting singularities are not of the observable type. 

The first threshold is located at M 2 , which is positive and larger than the free quark 
mass squared m 2 . For massless quarks, Mi remains positive. 

To analyze the solution in x-space, we define the Fourier transform of the decomposi- 
tion (13~TD : 

S(x) = H.dFx(r) + F (r), r = (3.6) 
We introduce modified functions Fx(r) and F (r) such that 

S(x) = ±-(P^Fx{r) + F (r)), F^r) = -2n^l, F (r) = 2nF (r). (3.7) 
ztx v r > or 

From the Fourier transforms of Eqs. (13.21) and ( 13.31) . one obtains: 



= o- E ^ e " M " r ' F oW = — £ (-l)" +1 6 n e~ M "'-. (3.8) 

2(T n=l 2(T n=l 

In terms of the functions Fx and F , Eq. (12. 3p decomposes into two equations (here 
written for r 7^ 0): 

3 ~ err 3 f /"^ r ~ ~ ~ ~ n 

— (rF!(r)) + mrF (r) = — | J dX A 2 - A)r)iq(Ar) - F ((l - A)r)F (Ar) 

- d£ - l)r)A(^) + F ((e - l)r)F„(£r)] }, 

(3.9) 
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— (F (r))+mF 1 (r) = — | jf dX A 2 [^((1 - A)r)F (Ar) - F ((l - A)r)iq(Ar) 



(3.10) 



Replacing in them Fx and _F with their expressions (13. 8p and considering each exponential 
function with a given M n as an independent function for r ^ 0, one obtains the following 
three sets of consistency conditions, valid for all n — 1, . . . , oo, n not summed: 

M„ = (-l)- m + !£(!- (-ir«) IJ ^A_ (3.U) 

1 = -£ (1 - ( - 1)M,) (ii4^?" (3 ' 12) 
£ (1 - ( - ir+ * ) W A ^) + ? (1 + (M^) = ° (3 ' 13) 

Equations (13. lip and ( 13 . 1 3[) are consequences of both of Eqs. (13.91) and (13.101) . while 
Eq. (13.121) results from Eq. (13. 9p only; its components play a particular role in the 
limit r — > 0. In that limit, the right-hand side of Eq. ( 13. 12ft . which is inserted in 
the series of the exponential functions, gives a vanishing contribution for antisymmetry 
reasons, while its left-hand side contributes globally as a divergent constant reproducing 
the delta-function of the right-hand side of Eq. (12. 3p . Equations (13 . 1 1 [) - f |3 . 1 3[) determine 
completely the mass parameters M n and m n , up to a global normalization factor fixed 
by the inhomogeneous part of Eq. ( 12. 3p . Details of the resolution method are presented 
in the appendix. Equations ( 13. lip and ( 13.121) mainly determine the parameters M. The 
parameters fn are asymptotically negligible in front of the Ms and they play a secondary 
role in the above equations. They are mainly determined through Eqs. (I3.13p . From the 
latter equations and from Eqs. ( 13.111) . one derives the following ones: 

which are more appropriate for an iterative determination of the Wis. 

Finally, we study the effect of the inhomogeneous part of Eq. (I2.3p . From Eqs. ( 13. 9p 
and ( 13.1 Op one finds that the interaction dependent parts, represented by the right-hand 
sides, vanish in the limit r = 0. In momentum space this means that the asymptotic 
behavior of S(p) is governed, for spacelike momenta, by the free propagator, as can be 
deduced from the inhomogeneous part of Eq. ( 12.31) . The asymptotic behavior of the scalar 
function F Q , however, crucially depends on whether the quark is massless or not [13] . 

We first check the behavior of Fi from Eq. (13.21) . Taking \p 2 \ large and neglecting 
in the right-hand side the M 2 s in front of p 2 , one finds a divergent coefficient which is 
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Friv 2 ) .= -, (3.15) 



an indication that the asymptotic behavior of F\ is slower than (— p 2 )~ 3 ^ 2 . To evaluate 
the asymptotic tail of the series, one can replace b n by its asymptotic behavior a 2 /M n ~ 
o 2 1 \foTm and convert the series into an integral. One then obtains for the asymptotic 
behavior 

i 

c p 2 

in agreement with what is expected from the inhomogeneous part of Eq. (12.31) . In x-space, 
Fi(r) behaves for spacelike x and small r as (l/(4n)) ln(l/(cxr 2 )) [16] and i*i(r) as 1/r. 

For the study of the asymptotic behavior of F , for m ^ 0, one again sticks to the 
asymptotic tail of the series which provides the dominant behavior. For this, one com- 
bines in the corresponding alternating series two successive terms, which leads to the 
factorisation of the parameter m. Proceeding then as for F±, one finds, as expected, 

F (p 2 ) = % m^O. (3.16) 

p 2 — >— oo P 

In x-space, F (r) behaves for spacelike x and small r as (m/(47r)) ln(l/(ar 2 )). 

When m — 0, the mass dependent terms in the left-hand side of Eqs. ( 13.91) and (13. 101) 
disappear. Eq. (13. 10p implies that J^(F (r)) vanishes at r = 0, which leads to a constraint 
on the masses and residues: 

oo 

£ (-l) n M4 = 0, m = 0. (3.17) 

n=l 

-F (0) is then a constant. Analyzing the behavior of the right hand-side at small r, one 
finds for F (r) the following behavior: F (r) = F (0) 1 — (ar 2 / (An)) ln(l/(crr 2 )) +0(r 2 ). 
By Fourier transformation to momentum space, the first term produces a delta-function 
that does not contribute to the asymptotic behavior; it is the logarithmic term that 
provides the dominant term in the asymptotic region [16] (r 2 = — x 2 , Eq. ( 13. 6ft ): 

F«(P 2 ) = - M " F ?*r 0) = %M, ™ = 0, (3.18) 

where in the second equation we have introduced the one-flavor quark condensate 
(?(0)V(0)> = -2iV c F (x = 0). 

We have thus shown that expressions ( 13.ip -( l3~3l) . or equivalently (I3.6p - fl3.8p . are solu- 
tions of Eq. ( 12.31) . provided the constraints ( 13 . 1 1 p - ( T37T3T) are satisfied. 

On practical grounds, there remains to solve the latter equations for the parameters 
M n and m n (n = 1,2,...). The properties of those equations are studied in the appendix, 
where also details about their numerical resolution are presented. We have calculated, for 
several values of the quark mass, m = 0., 0.1, 1., 5., 20., in mass unit of ycr/ft — 240 
MeV, the first 900 values of the parameters M n and rn n , with better accuracy than 10~ 4 
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for the first 400 values. For massless or light quarks, the asymptotic values (I3.5p are 
reached within a few per cent accuracy after a few dozens of ns. With increasing quark 
masses, the asymptotic limit is reached more slowly. Approximate formulas for the next- 
to-leading terms in the asymptotic behavior of M n are presented in the appendix. We 
display in Tables [T] and H] a few typical values of M and m. 



n 


1 


2 


3 


4 


5 


6 


225 


226 


m = 


1.161 


3.043 


4.301 


5.287 


6.127 


6.868 


46.96 


47.06 


m = 0.1 


1.248 


3.103 


4.344 


5.323 


6.158 


6.897 


46.96 


47.06 


m = 1 


2.053 


3.706 


4.823 


5.745 


6.536 


7.246 


47.03 


47.14 


m = 5 


5.833 


7.059 


7.898 


8.642 


9.287 


9.887 


47.92 


48.02 


m = 20 


20.59 


21.41 


21.97 


22.48 


22.93 


23.36 


55.61 


55.70 



Table 1: A few values of M, in mass unit of ya/n, for several values of the quark mass 
m. 



n 


1 


2 


3 


4 


5 


6 


225 


226 


m = 


1.052 


0.217 


0.282 


0.147 


0.191 


0.120 


0.022 


0.021 


m = 0.1 


0.969 


0.207 


0.282 


0.143 


0.192 


0.117 


0.022 


0.021 


m = 1 


0.553 


0.145 


0.257 


0.110 


0.189 


0.094 


0.023 


0.019 


m = 5 


0.179 


0.056 


0.143 


0.051 


0.125 


0.048 


0.025 


0.016 


m = 20 


0.049 


0.016 


0.047 


0.016 


0.046 


0.016 


0.022 


0.010 



Table 2: A few values of m, in mass unit of ya/vr, for several values of the quark mass 
m. 

We present in Fig. Q] the function iF Q for spacelike p and in Fig. [2] its real part for 
timelike p, and similarly for the function iF\ in Figs. [3]and|U for m = 0. Their imaginary 
parts (the spectral functions multiplied with 7r) have behaviors that resemble those of 
the real parts in the timelike region, except that they are null below Mi and that the 
infinite singularities begin just from the right of the Ms. A relationship between iFi(0) 
and ^-Fo(0) exists; it is presented in the appendix. 

The x-space functions F Q (r) and F±(r) [Eqs. (I3.7p -( l3~8]) ] are positive and monotone de- 
creasing functions of r for spacelike x. For timelike x, the two functions are exponentially 
oscillating. 
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iF 




Figure 1: The function iFq for spacelike p, in mass unit of J a /it, for m = 




iF 1 




Figure 3: The function iF± for spacelike p, in mass unit of J a /it, for m = 



fte(iFi) 



1 H 




M| 



Mr 



M 7 2 



Figure 4: The real part of the function iF\ for timelike p, in mass unit of ya/n, for 
m = 0. 
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4 Momentum space analysis: the instantaneous limit 



A direct resolution of Eq. (12. 3p in momentum space is complicated because of the presence 
of the additional integrations upon the variables A and £. Nevertheless, a two-step reso- 
lution, using the analyticity properties resulting from the presence of spectral functions, 
is possible. This section and the following one describe the corresponding procedure. 

The starting point is based on the observation that Eq. (12.31) is quasi-local in x, in the 
sense that in the right-hand side of the equation x undergoes only dilatation operations. 
Therefore, projections of the Green's function on particular subspaces still give rise to the 
same equation (for x ^ 0). Of particular interest is the instantaneous limit x° = 0. In 
that case, the variable r of Eq. (I3.6P simply becomes the modulus of the space component 
x 1 : r = \x l \. The functions Fi(r) and F (r), introduced in Eqs. (13.71) . keep the same 
definitions and Eqs. (13. 9p and (I3.10p remain unchanged in form with the only modification 
that r has become a one-dimensional variable. 

We pass now to momentum space with one-dimensional Fourier transformation. The 
momentum variable k will designate the modulus of the space component p\ of the two- 
dimensional vector p: k = \pi\. The Fourier transform of Fq(v) will be defined with the 
cosine function, while that of F\(r), which is defined as a derivative of F\(r) with respect 
to r [Eq. (13.71) ]. with the sine function. We have the definitions 

/•oo roo 

f (k)=2 dr cos(kr)F (r), fi(k) = 2 1 dr sin(kr) F x (r), r = \x l \, (4.1) 
Jo Jo 

„ (1A [°° dp f (k) r°° dp U(k) 

(4.2) 

Replacement of -F\(r) and F (r) in Eqs. (13 .9p and (I3.10p by their one-dimensional 
Fourier transforms allows the integration with respect to the variables A and £. One ends 
up with the equivalent equations for f\{k) and fo{k): 



u dfi(k) df (k) 
k — — h m- 



dk dk 2w 2 Jo 



~ UM ^k~ ~ fo{q) ^k~ ) J^TW 
+(/i(9)-air + /oC?)-^-) j^Tky } - -Ji 

(4.3) 



kh{k)-mf x {k) = ^J™dq{(f l (q)f (k) + fo(q)h(k))-^ 



<fi(q)fo(k) ~ fo{q)h{k)) ^— ^ }. (4.4) 
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Equations (14 .3p and (I4.4p display singular kernels. Their effects on the solutions can 
be studied by expanding the variable q around k in the functions fi(q) and fo(q). It turns 
out that Eq. (14.41) is finite, while Eq. (14.31) is in general divergent. Finite solutions exist, 
provided that they satisfy the constraint fi+fo =constant. The value of the constant can 
be determined from the inhomogeneous part of Eq. (I2.3P and the property that f\ and 
fo should behave asymptotically as in a free field theory, i.e., f\{k) — > n and fo(k) — > 
for k — > oo. This fixes the constant to n 2 : 

/?(*) + /o(*) = ^ (4.5) 
Parametrizing the functions j\ and fo by means of an angle (p(k) as 

fi(k) = n cos (p(k), f (k) = 7rsm<f(k), (4.6) 

Eq. (14.31) reproduces Eq. (14.41) in addition to its last term (the derivative of the delta- 
function). The compatibility of the two equations then requires the vanishing of the 
coefficient of the latter term, i.e., of /i(0): 

/i(0)=0, or, ^(0) = |. (4.7) 

The final equation is then: 

a r f°° 1 1 

ksinp(k) -mcosp(k) = — [ cos^/c)^ dq sin <p(q) ( + + - _ ) 



f°° ll-i 
+ Bm<p(k)J Q dqcosipmj^-j—w)}- 



{q + k) 2 (q-kf 

(4.8) 

Let us emphasize at this point that the boundary condition ( 14 .7p . which is a necessary 
condition for the existence of solutions to Eqs. (I4.3p - (l4.4p . eliminates the free case ip = 
when m = from the set of possibilities, which otherwise would fit Eq. (14. 8p . This 
conclusion can also be reached directly in x-space, by inspecting Eq. ( 12.31) . 

Equation (14. 8p is identical to the one satisfied by the ordinary quark propagator in 
its instantaneous limit in the axial gauge [TTJ, [HI US]- Therefore, the 2PGIQGF and the 
ordinary quark propagator in the axial gauge coincide in their instantaneous limits. 

More generally, Eq. (14. 8p is analogous to the self-energy equation obtained in the 
instantaneous limit in the four- dimensional theory, using Coulomb gluons for the descrip- 
tion of the confinement mechanism [20l [2TI 1221 1231 [24] . Here, however, as well as in the 
axial gauge of the two-dimensional theory, the ordinary quark propagator is a divergent 
quantity due to the infra-red singularity of the gluon-propagator, and is made finite with 
the aid of a regularization procedure. It turns out that its instantaneous limit is free of 
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divergences and satisfies a finite equation [22J. This is the reason why the interface be- 
tween the 2PGIQGF, which in general is expected to be an infra-red finite quantity, and 
the ordinary quark propagator in a non-covariant gauge, exists only in the instantaneous 
limit. 

From Eq. (14 .8|) . one deduces that asymptotically, as k — > oo, sirup behaves as m/k 
when m^fl and as If ^ sinyj(g) when m = 0, in agreement with the results obtained 
in Eqs. ( I3TT6]) and f|37T8]) . 

Equation ( 14.8}) can be solved with the same methods as its four-dimensional analogue 
[221 125]. The solution for sin</>, for massless quarks, is presented in Fig. |5j 

1.25 i 1 1 1 1 1 1 



sin if 




1 1 1 1 1 — 1 1 

1 2 3 4 5 6 

k 

Figure 5: The function sin Lf, in mass unit of Ja/ir, for m = 0. 



5 Analyticity properties in the instantaneous limit 

The numerical knowledge of the instantaneous limit of the functions Fi(p) and F Q (p) does 
not allow an immediate reconstruction of the whole covariant functions and in particular 
of their spectral functions. Nevertheless, the existence of underlying analyticity properties 
allows us to reach that aim by means of an analytic continuation. This section is devoted 
to that procedure. 

To simplify the presentation, we shall assume that the series of spectral functions 
(11.21) sum up, by means of integrations by parts, into the lowest-order denominator. If 
the spectral functions possess strong singularities that forbid the above operations, one 
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can bring outside the integrals the appropriate number of derivation operators. The 
spectral representation of the functions F\ and F then becomes: 



CO 



W) = W ds 2 P1(S ]. F (p>) = i nds . 2 Po{s ]. v (5.1) 

JO (p 2 - S + IE) Jo (p 2 - S + IE) 

Taking the instantaneous limit of the latter functions [Eqs. ( 14.2ft and ( 14.61) ] . one finds 

cosy(fc) = k f°° ds , sin^(fc) = r ds . (5.2) 

Jo \Jk 2 + s Jo \Jk 2 + s 

In principle, the above integral equations could be inverted with the successive uses 
of the Hankel transform of order zero and the inverse Laplace transform [26]. However, 
the latter transform necessitates knowledge of analytic expressions. Also, the direct cal- 
culation of the product of the two transforms gives a singular expression, since the Bessel 
function of order zero is an analytic function. These difficulties can, however, be circum- 
vented with the use of an analytic continuation method. 

The key observation is that because of the properties of the spectral functions, assumed 
to be concentrated on the positive real axis, the instantaneous limit functions cos<£> and 
sin ip are themselves analytic functions of k in its whole complex plane with possible sin- 
gularities located on the imaginary axis. We are then entitled to continue those functions 
to complex values of k and in particular to the imaginary axis. 

Let us consider for definiteness the function sin ip. One can express it as the Laplace 
transform of some other function, which will be denoted 

sincp(k) = [°° dyh (y)e-y k . (5.3) 



The exponential function e yk can be expressed as the Hankel transform of order zero of 
the function l/y/k 2 + s: 

e~y k = y - rd8-=±=J (yy/s), (5.4) 
2 Jo y/k 2 + s 

where Jo is the Bessel function of order zero. Replacing in Eq. (15.31) the exponential 
function by its expression (15.4j) . one arrives at the identification 

1 r°° 

Po{s) = ~ dyh (y)yj (yy/s). (5.5) 



Finally, using for the Bessel function the representation Jq(z) = (l/7r) J^i dte tzt j\f \ — t 2 
and continuing Eq. ( 15. 3 j) to complex values of k, one obtains the formulas 



1 /V- , Id 

Po{s) = 



/•V s 1 (J/ 

\ du—^== — [sin <p(-iu) -sin <p(iu) 

2t7T Jo y/S — U 2 OU K 

^- [ duy/s — u 2 ^—( sintpi— iu) — simpiiuj). (5.6) 

2in os Jo ou K J 
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Assuming that the real part of sin p is an even function on the imaginary axis and 
the imaginary part an odd function, Eq. ( 15. 6 j) means that the spectral function p is 
determined from the imaginary part of the function sin ip on the imaginary axis. 

A similar calculation as above can be repeated for the function cos tp. One obtains 

1 fV* 1 d rl i 

Pi{s) = — du = - (cosip(-iu) + coscp(iu)) 

Ztx Jo a/s — u ou l u J 

2 d rV^ d r l i 

= — — / duy s — u 2 — —(cos ip(—iu) + cos <p(iu)) . (5.7) 
2n os Jo ou l u J 

The spectral function p\ is then determined from the real part of cos p on the imaginary 
axis. 

To have an elementary check of the consistency of Eqs. (15. 6p and (15. 7p . let us assume 
that the spectral functions are null for values of s below a positive mass-squared M 2 . 
Then, for values of k smaller than M, one has from Eqs. (15.21) the expansions sinp(k) = 
1 + ak 2 + 0(k 4 ) and cosip(k) = bk + 0(k 3 ). The analytic continution to the imaginary axis 
gives: sin p(±ik) = 1 — ok 2 + 0(k A ), cos <p(±ik) = ±ibk + 0(k 3 ), from which one verifies, 
through Eqs. (15. 6p and (15 .7p . that p and pi, as expected, are null in the considered 
domain. 

The evaluation of expressions ( 15. 6ft and ( 15.71) necessitates the analytic continuation 
of the functions sinp and cosp to the imaginary axis of the complex variable k. This 
can be achieved by using the integral equation (14. 8p . in which the integrands sinp(q) 
and cos p(q), which lie on the positive real axis, will now be assumed numerically known 
from the resolution of the equation for real positive values of k. The validity of that 
continuation is justified by the fact that Eq. (14. 8 p is the instantaneous limit of covariant 
scalar equations relative to the functions F\ and Fq. The singularities of the integrals arise 
from the terms l/(q — k) 2 , which actually contribute, by combining the various terms, as 
a principal value integrand. (The corresponding definition results from the Fourier sine 
or cosine transforms of the terms of the one-dimensional limit of Eq. (12.31) .) The latter 
has a well-defined meaning as a generalized function [16]. One has: 

[P designates the principal value.] The first term in the right-hand side of Eq. (15.81) is 
the boundary value of an analytic function and can readily be continued to the complex 
plane. According to the sign in front of ie, the continuation will be done to the lower 
half-plane (minus sign) or to the upper half-plane (plus sign). The delta- function can be 
integrated with respect to q (actually, preceded by an integration by parts) and yields a 
term proportional to -^p(k), which is an analytic function of k and can be continued to 
the complex plane. 
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We shall mainly be interested by the limit of k to the imaginary axis (with a small 
positive real part) and therefore shall directly consider this case: 



k^r]±iki, r] > 0, r] ~ 0, ki > 0. (5.9) 

[For simplicity, rj will henceforth be omitted from the formulas.] With the above defini- 
tions, the continuation of Eq. (14. 8p to the imaginary axis takes the form 

(j d ex 
±iki sm(p(±iki) — m cos ip(±iki) H — — — <f(±iki) ^fi—<f'(0) sinip(±iki) 

2 ut\"i 2 

= cos <f(±iki) g ± (ki) =F isin<f(±iki) g 2 (ki), (5.10) 
where <f'(0) = 9ip d ^ \k=o, and g\ and g 2 are the following integrals: 

U roc (q2 

9l (h) = -J o ^7^-^sin^g), (5.11) 

g 2 {h) = - <ig arctan(— ) — -cos^(<?) 
it Jo q oq z 

U (j r°° 1 d 

= -<p{0) + -k i J Q dqj^^^cos^q). (5.12) 

The first expression of g 2 is useful for small values of ki, while the second one for large 
values of k^ one has #2(00) = f y?'(0). g\ and g 2 are in general negative. In magnitude, |gi| 
is a monotone decreasing function from a finite value at ki — down to zero at infinity; 
|<?2| is a monotone increasing function from zero at ki = up to a\if'(0)\/2 at infinity 
(f'(0) is negative). 

The functions gi(ki) and g 2 (ki) being known from the resolution of Eq. (14.81) . we find 
that the continuation of the latter equation to the imaginary axis has yielded a first-order 
differential equation in if, which can be analyzed in turn. 

Decomposing ip into a real and an imaginary part, f{k) = f r (k) + ifi(k), Eqs. (15.1 Op 
decompose into two real equations: 



±fcj sin (p r (±iki) cosh ifi(±iki) + m sin if r (±iki) sinh ifi(±iki) 
+-T^(Pi(±iki) =F -<f' (0) sin ip r (±iki) cosh. cpi(±iki) 
= — sin if r (±iki) sinh (fi(±iki) gi(ki) =F sin ip r (±iki) cosh ipi(±iki) g 2 (ki), 

(5.13) 

^k{ cos if r (±iki) sinh (fi(±iki) — m cos <f r (±iki) cosh ifi(±iki) 

(7 d (j 
+-—ipr(±ih) ± -<f'(0) cos f r (±iki) sinh ifi(±iki) 

= cos f r (±iki) coshifi(±iki) g±(ki) ± cosf r (±iki) sinh ifi(±iki) g 2 {ki). 

(5.14) 
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The boundary values to be imposed are y? r (0) = n/2 and <Pi(Q) = 0, i.e., (p(0) = n/2 [Eq. 

The above equations can be solved starting from a domain containing the origin, to 
implement the boundary conditions. It is found that the singularities of the functions 
sin<^ and cosy? on the imaginary axis are represented by an infinite set of simple poles 
with locations at k — ±iM n , n = 1,2, ... . Considering for instance the positive imaginary 
axis, one can divide it into domains, each one lying between two successive poles, the first 
domain lying between the origin and the first pole. In the first domain, one has ip r = 7r/2; 
in the second one, <p r = — vr/2; in the third one, ip r = tt /2 (modulo 2ir); ip r thus alternates 
between the constant values ±7r/2 when passing from one domain to the other. 

In each domain, ipi is a monotone increasing or decreasing function of ki (taking into 
account the signs and values of the functions g± and g 2 ). In the first domain, the solution 
for ifi is given by the equation 

crarctan(tanh( V9 '^ )) = J - k[ + ^<p'(0) - g 2 (k[) - (m + gi{k[)) tanh^fc-)) , 

(5.15) 

which can be solved by iteration. [We henceforth omit in the argument of p>i the factor i.\ 
The position of the singularity M\ corresponds to that value of ki for which the right-hand 
side reaches the value air/ 4; at that position, tpi becomes logarithmically infinite. One 
has then the equation 

I dk' l (k' l --p'(0) + g 2 (k' 1 ) + (m + gi (k' i )) tanh^(fcO) = a-. (5.16) 

In the (n + l)th domain {n > 1), M n < ki < M n+ i, the corresponding equations are: 

aarctan(tanh(^^)) = (-l) n {a- 

+ [ ki dk>( -k> + V(0) - M + 9i{K)) tanh^(fcO)}, 

(5.17) 

Jm d K{K-^\0)+g 2 (k^ + (m + g 1 (k^)tanhip l (k' i ))=a-. (5.18) 

The residues of sin ip and cos ip at the position of their poles are found by expanding 
in the various expressions ki around M n . One finds: 

sm<p(±iki) ~ (-l) (n+1) -7TT —, : TTTT \ r, (5.19) 

yi fc^M„ ; 2{M n + m n + {-l) n m)(M n -ki±ieY 1 ; 

1 1 

cos(p(±iki) ~ ± i — — — ; r tttt ; r, (5.20) 

^ k^M n 2{M n + m n + {-iym){M n -k l ±iey 1 ; 

m n = ~tp\0) + g 2 (M n ) + (-l) n gi (M n ). (5.21) 
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The real parts of sin tp and i cos (p for k on the positive imaginary axis are shown in Figs. 
[6] and [71 They are even and odd functions of ki, respectively. 




The spectral functions po and p\ are obtained from Eqs. (15.6p and ( 15.7ft . One finds: 
PO W - 2 2^ lj + _ + ( _ 1)Bm) (a _ M 2 )3/2 , (5.22) 

Pl[) - 2&{M n +m n + (-l)»m){8-M*)*/*- { - ] 

The functions siny? and cosy? on the positive real axis are reconstituted from Eqs. 
05 .2p . [Since the ps have strong singularities at the various thresholds, one has to use for 
them the first version of equations (15.61) and (15.71) and make an interchange in the order 
of integrations between u and s.\ One obtains: 



oo 



sinp{k) = -a E(- 1 ) n 777 ^ > ^ 7 775 ^ (5-24) 

00 1 1 

cosp(k) = ak V* — — ; — — -— — — — . (5.25) 

~i (M n + m n + (-l)»m) (fc 2 + Ml) v ; 

As a check, one verifies that the curve of sirnp(k) obtained from Eq. (15.24ft with the 
numerical values of the mass parameters M n and m n coincides with that obtained from 
the direct resolution of Eq. (I4.8p [Fig. [5] . 

One can go further, by transforming the informations coming from the instantaneous 
limit into an infinite set algebraic equations. First, going back to the defining equations 
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Figure 7: The real part of the function i cos tp on the positive imaginary axis. 



of the parameters m n [Eqs. (I5.2ip ] and calculating now the functions g\ and <?2 with the 
use of the expressions (I5.24p and (I5.25P of sin tp and cos <p on the real axis, one finds the 
set of equations (I3.14p allowing the calculation of the ms for known Ms. 

Second, using expressions (I5.24p and (I5.25P in Eq. (14. 8p . one finds, as consistency 
conditions Eqs. (13. lip and (I3.13p . Equation (13 . 12[) is implicit here, since odd-index and 
even-index parameters succeed to each other. [Their definition is not a mere artifact, 
their defining equations being different, cf. Eqs. (I5.17P and (I5.2ip .] The equation sin 2 <p + 
cos 2 (p = 1, which is a non-trivial constraint for the expressions (I5.24p -f l5.25p . leads also to 
the previous consistency relations. One thus finds at the end the same algebraic equations 
as those found in the x-space analysis. 

Finally, the expressions (I5.23P and (I5.22p of the spectral functions allow us to recon- 
struct the completely covariant functions F\ and F , given in Eqs. (13. 2p and (13. 3p . 

This completes the study and determination of the quark Green's function in momen- 
tum space. 



6 The quark condensate 

A quantity of interest is the quark condensate, which, in the one-flavor case and for 
massless quarks, has the expression [Eq. (13. 6p ] 

(^(0)^(0)) = -ti C:Sp S(m = 0,x)\ x=q = -2N c F (m = 0,x = 0), (6.1) 
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tr c tSp meaning the trace operation taken on color and spinor indices. This quantity, being 
the vacuum average of a local gauge invariant operator, is also calculable from the ordinary 
quark propagator and does not need for its evaluation the recourse to gauge invariant 
Green's functions. Nevertheless, it is interesting to compare, as a check, its value obtained 
in the present approach to that obtained in previous approaches. 

The quark condensate is calculated from the series expression of the scalar function 
F (x) [Eqs. (j3T?p . (|53jl and lETip ] for m = 0: 

(^(0)^(0)) = -f- £ (-1)" — ~ -293 MeV, (6.2) 

the numerical value being obtained with iV c = 3 and yfa = 425 MeV. [The asymptotic 
tail of the series is evaluated using the asymptotic expression of M n and neglecting the 
ms.] 

The above numerical value is in agreement with that obtained in the axial gauge from 
the integral of the function siny? for massless quarks [18] ({ipip) — —295 MeV) and with 
that obtained in the light-cone gauge using operator product expansion methods [23 |2B] 
= -N c ^Ja/(6n) ~ -294 MeV). 



7 Conclusion 

The study of the gauge invariant quark Green's function in two-dimensional QCD in 
the large- N c limit by means of the integrodifferential equation it satisfies has allowed us 
to determine the properties of its spectral functions. Those functions directly probe the 
contributions of quarks and gluons in a sector where no color-singlet states can contribute. 
The main results can be summarized as follows. 

1) The spectral functions are infra-red finite and lie on the positive real axis of the 
variable p 2 in Minkowski space (p timelike). No singularities in the complex plane or on 
the negative real axis have been found. This feature underlines the fact that quarks and 
gluons contribute to the spectral functions like physical particles with positive energies. 

2) The singularities of the spectral functions are represented by an infinite number of 
threshold singularities, characterized by positive real masses M n (n — 1, 2, . . .) extending 
from a lowest positive value up to infinity. The corresponding singularities near the 
thresholds are of the type (p 2 — M 2 )~ 3 ^ 2 and therefore are stronger than simple pole 
singularities, a feature which reflects the difficulty in observation of quarks as free particles. 

3) The threshold masses M n remain positive in the limit of massless quarks and main- 
tain the scalar part of the Green's function at a non-zero value. 

The present method of investigation can also be applied to four-dimensional QCD, 
in which case, however, several difficulties arise. First, on practical grounds, the inte- 
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grodifferential equation of the gauge invariant Green's function can no longer be solved 
in an exact way, because of the presence of an infinite number of kernels and of the 
different contributions of short- and large-distance interactions. Second, the presence of 
short-distance interactions induces divergences which necessitate an appropriate treat- 
ment through renormalization. Third, the Wilson loop averages are no longer known 
exactly. 

Another domain of investigation is the four-point gauge invariant Green's function, 
which is naturally related to the bound state equation of mesons. 
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A Determination of the spectral function parameters 

This appendix is devoted to the determination of the threshold masses and factors that 
appear in the spectral functions. In the first part, we study the mathematical properties 
of the corresponding equations; in the second part, the asymptotic behaviors of the mass 
parameters M n and m n for large n are found and in the third part, details of the numerical 
resolution of the equations are given. 

A.l Mathematical properties 

We notice that Eqs. (13. lip and (I3.12p separate into two complementary sets according to 
the parity of the index n: 




n = 0,1,2, ... 




oo 



1 





(A.2) 



oo 



M> 




n = 0,1,2,... , 



(A.3) 



oo 



b 




Vt-i 




(A.4) 
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The odd-index Ms are thus determined by equations where the even-index Ms play the 
role of background data. The roles of odd-index and even-index Ms are interchanged in 
the complementary sets of equations and m replaced by — m. 

To analyze the constraints imposed by Eqs. (1A.ll) -( lA~4~l) . we may begin by considering 
Eqs. (lA.ip as representing the intersection condition of two curves. Let us define the two 
functions 

oo ^ 

M x ) = H (rr m Y» » f+2(x)=x-m, (A.5) 

q =0 \ X - M 2g+2) 

for x real and where the even-index Ms are assumed to be known. Then, the solutions 
of Eqs. flA.ip represent the intersection points of the straight line f + 2 with the curve 
f + i. We further notice that the function / +1 has an infinite number of double poles at 
the even-index masses with positive residues. Any solution M2 n +i should lie between 
two successive even-index poles. However, in general, the equality f + 2 = f+i might 
not lead to such a solution. First, the straight line f + 2 might not intersect at all the 
curve in which case no solution would exist. Or, second, it might cut f + \ at two 
different points between two successive even-index poles. In that case, we would meet 
instability, since, considering then the complementary sets of equations with the roles of 
odd and even indices interchanged, and with the number of odd-index solutions doubled, 
we would generate new solutions of even-index Ms (doubled in number), and so forth. 
We therefore guess that the only stable and acceptable solutions to Eqs. (1A.1I) are those 
corresponding to the situation where the straight line f + 2 is tangent to the curve f + \ 
between two successive even-index poles. This is precisely the content of Eqs. (1A.2I) . 
Written as functions of x, they take the form 

df+2(x) = df +1 (x) ^ A g , 

dx dx ' 

indicating that the straight line / +2 should be tangent to the curve f + i at the positions of 
the solutions (cf. Fig. [8]). Of course, this condition could not be satisfied if the even- index 
Ms and ms were arbitrary. The even-index Ms should themselves be solutions of Eqs. 
(1A.3|) and (IA.4j) ensuring the final consistency of all equations. 

While the geometrical interpretation of Eqs. (lA.lj) - (lA.4j) is very suggestive of the 
existence of acceptable solutions, one might still question the mutual compatibility of the 
latter equations. Assuming that the mass parameters m, or equivalently the residues b, 
are mainly determined from Eqs. ( 13.1 3D or ( 13. 14ft . one notices that the system of equations 
(1A.1|) - (1A.4|) represent four sets of equations for two sets of unknowns, M 2n+ i and M 2n+ 2, 
n = 0, 1, 2, . . . . It is then necessary to check the compatibility of the above equations. 

To that end we define two meromorphic functions h+(z) and h-(z): 

OO L 

°2q+2 

„ (z - M 2 , +2 ) 
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Figure 8: Graphic representation of Eqs. (1A.1I) and (1A.2I) for m — 0. The solutions, 
for odd-index parameters, for given even-index parameter solutions, correspond to the 
contact points of the straight line with the curves. 



=0 V Z ~ M 2q+l) 



Or 



The singularities of h + are double poles located at the even-index masses, while those 
of h- are double poles located at the odd-index masses. The residues b n are defined 
in Eqs. (13.41) and are assumed to be finite quantities. Equations flA.lj) take the form 
h + (M2n+i) = 0, n = 0, 1, 2, ... , which means that h + has zeros at the positions of the 
poles of h_. Equations (I A. 21) take the form dh +^ = Q which means that the 

H " " dz z=M 2n+1 

zeros of h+ are double. Therefore, the product h+h- is finite at the poles of A similar 
argument, using now Eqs. (1A.3I) and (1A.4I) . also holds at the poles of h + . In conclusion, 
the product h + (z)h_(z) is an analytic function in the whole complex z-plane. According 
to Liouville's theorem [29], it is a polynomial of second degree: 

h+(z)h- (z) = z 2 -m 2 + — , (A.9) 

7T 

where the right-hand side has been determined from the asymptotic behaviors of h + and 
h- . It is to be emphasized that the formal calculation of the product of the series flA.7j) 



and flA.8j) by exchanging the orders of summation in the double sum would miss the last 
constant term in the right-hand side of Eq. flA.9j) ; this is due to the weak convergence of 
the individual series. 
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Equation flA.9[) indicates that if h + is known, then h_ is completely determined from 
it. It thus provides the proof of compatibility of Eqs. (1A.1|) - (1A.4|) . For suppose that 
Eqs. (lA.ip and (1A.2|) are considered on formal grounds as two independent sets of equa- 
tions determining the two sets of unknowns M2 n +i and M2 n +2- Then the two remaining 
equations (1A.3j) and (1A.4j) would not bring any new information or constraint since their 
content is already summarized in Eq. flA.9j) . On practical grounds, it is still preferable 
to determine by means of iterative procedures the odd-index masses from Eqs. (1A.1I) and 
flA~2l and the even-index masses from Eqs. (1A.3|) and (IA.4I) . 



Relations, involving the residues b n , are obtained by taking in Eq. ( 1A.9I) the limit of 



z to one of the pole positions of h+ or For odd-index poles, for instance, one finds: 

oo L n 

3^.E , M W M ■' + f. (A.10) 

q=0 [ M 2n+1 — M2q+2) ^ 

with a similar relation for the even-index poles. The procedure can be continued with 
higher-order derivatives at the poles, leading to new relations for the b n s, but it is not 
evident whether the derivation operation inside the series remains still valid at higher 
orders. Equations ( 1A.9I) and ( 1A.10I) can be used as consistency checks for the solutions 
that are found. 

We next establish Eq. (I3.14p . The derivation follows similar lines as above. We define 
two meromorphic functions f+(z) and f-(z): 

/ + W - i + £ ) ■ < A - n) 

q=0 \ Z - M 2q+2) q=0 \Z + IVhq+l) 

/-(*) - t TT^Sh-, + t TT^T^-v (A- 12) 



q=0 \ Z ±y± -lq+l) q= 

The function / + has poles on the positive real axis at even-index masses and on the 
negative real axis at minus the odd-index masses. /_ has analogous properties with the 
roles of even-index and odd-index masses interchanged. / + and /_ satisfy the following 
relation: 

/_(*) = -/+(-*). (A.13) 
Equations f 1 3 . 1 3 j) can be rewritten in the form 

/+(M 2n+1 )=0, /_(M 2n+2 ) = 0, n = 0,1,2,... . (A.14) 

This means that / + has zeros at the poles of /_ and /_ has zeros at the poles of / + . 
Therefore, the product f + (z)f_(z) is free of singularities and is analytic in the whole 
complex z-plane. According to Liouville's theorem [29], it is a constant, which could be 
determined from the asymptotic behaviors of /+ and /_ at infinity: 

f + (z) f_(z) = -a 2 . (A.15) 



M 2q+ i) {z + M 2q+2 \ 
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[Because of the weak convergence of the series defining the functions f±(z), it is not pos- 
sible to obtain the value of the constant by formal calculation of the product f + (z)f-(z), 
the latter operation being invalid under the exchange of the orders of summation in the 
double series. Asymptotically, the functions f±(z) tend to imaginary constants.] 

It is then possible to calculate the derivatives of / + and /_ at the position of their 
zeros. Thus: 



dz 



a 



*M* n +x\fl(z) dz 1 



a 



(A. 16) 

z=M 2 „+i z->A/ 2 „+i V /?(z) dz J &2ra+l 

Calculating the expression of the left-hand side and then using Eq. flA.lj) . we obtain the 
equation: 



>2n+l 



M 2n +i -m + Y^ 



J 2q+1 



q=Q 



2n+l 



which yields 



a 



J2n+1 



2n+l 



m 



!>■> 



9 + 1 



(A.17) 



(A.18) 



'q=0 {M 2n +i+M 2q+ i) 2 

Comparison with Eq. (13 .4p for odd n then provides the expressions 



m 2n +i 



E 



J 2q+1 



" Q (M 2n+1 + M 2q+1 ) 



(A.19) 



which correspond to the odd- index parts of Eqs. ( 13.141) . A similar derivation can also be 
done with the even-index parameters. 

Actually, the b n s could be defined with a different global multiplicative factor than 
in Eqs. (13.41) . The choice done there is the one which satisfies the correct normalization 
condition of the function Fi, as emphasized in Sec. |3j 

Coming back to the functions h±(z) [Eqs. ( 1A.7j) and (1A.8[) ] . it is possible to establish 
through them a relationship between the functions iF\{p 2 ) and iFo(p 2 ) [Eqs. (13. 1ft . ( 13. 2 h 
and (I3.3p ] at p 2 = 0. We also include for this analysis the series expansions of their 
instantaneous limits defined in Eqs. (TCI]), (02]), (IP)) . d4\7j) . ( l5T24"j) and ( 15T25|) . and which 
involve the function tp(k) parametrizing them (k is the modulus of the space component 
Pi of the vector p). Considering the particular values h + (0) and h_(0), the derivatives 
h' + (0) and h'_(0) and Eq. (1A.9j) . one finds the following relationships: 



iF (0) 



2m - 

mn 
a 

71 

27 " 



2a 



7T 



zF (0) 



4a 



r7T 



^(0) 



It 2 7T 

— (m — 
a a 



zF (0), 



1/2 



7T 



v'(o)[-^' 2 (o) 



TT 2 7T 

—lm — 
a a 



2) 



1/2 



(A.20) 
(A.21) 
(A.22) 



26 



where <f'(0) = , n . Elimination of <f'(0) yields a relationship between iFi(0) and 
iFo(O): 
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21 u v"/ 
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iFx(O) = — + - iF o (0) (^o(O)r - 2^F (0) + — . (A.23) 
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Numerically, one finds, for m = 0, in mass unit of Ja/n, from the series expansion 
-<V(0) = T,n=i K/Ml [Eq. (T5T25T) ] and the values of M n and m n , <p'(0) = -1.325, which 
predicts IFi(O) = 1.512 and iFq(0) = 1.528, in complete agreement with their direct 
calculation from the series (13 .2D and (13. 3L 

Other relations can also be obtained, like the expressions of the functions <p' and tp 
in the form of series and alternative expressions of the residues h n in the form of infinite 
products of factors. However, these results are not of direct practical relevance for the 
present calculations and will not be displayed. 



A. 2 Asymptotic expansions 

The asymptotic behaviors of the mass parameters M n and m n for large n can be studied 
from Eqs. ( 15. 18ft and (15.211) . If the quark mass m is itself large, then one has to distinguish 
two different regions: one for which aixn ^> m 2 and one for which awn m 2 . If m is 
small or null, then the first region may be identified with the region defined by the large 
n values. 

We first consider the first region (k 2 ~ ann ^> m 2 ). In the integral (I5.18p . tanh^j(fcj) 
varies from ±1 to =)=1 when ki varies between the two mass bounds M n and M n+ \. There- 
fore, it is the first term of the integrand that dominates and yields a recursion relation 
that can easily be solved leading to Eq. ( 13. 51) . [The other contributions of the integrand 
vanish at large ki] The corresponding behavior of m n [Eq. (13.51) ] is obtained from the 
contribution of the combination — §<£>'(0) + g 2 (A^n) to Eq. (15 . 2 1 D . g 2 being defined in Eq. 

ITED. 

To obtain an estimate of the next-to- leading term in the large n behavior of M n , it 
is necessary to specify in more detail the structure of the function tanh (fi(ki). Using 
the analytic continuation of the functions simp(k) [Eq. (I5.24p ] and cos</?(/c) [Eq. (I5.25D ] 
to the imaginary axis, it is easily seen that tanhyji(fcj) can be expressed in terms of the 
functions / + and /_ introduced above [Eqs. (1A.11|) and (1A.12j) ]: 

In order to evaluate the integral of mtanh <fi(kj) in Eq. (15 . 1 8j) . one can approximate the 
functions / + and /_ by the contributions of the pole terms in M n and M n+1 and take into 
account the contribution of the remainder, whose role is mainly to ensure the condition 
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that /+ (/_) vanishes at the positions of the poles of /_ (/+) [Eqs. flA.14j) ]. in the form 
of effective fcj-dependent residues at the above two poles. Thus, if we are integrating 
between M 2n and M 2n +i, we could use the approximations 

t b 2n (ki - M 2n+1 ) b 2n+1 (kj - M 2n ) 

J+Kl) (M 2n - M 2n+1 ){h - M 2n ) ' 1 tJ (M 2n+1 - M 2n )(h - M 2n+1 ) ' 1 • ; 

for M 2n < ki < M 2n+ i, and similar ones for the interval [M 2n+ iM 2n+2 \. 

With the above approximations, the calculation of the integral of tanh (fi(ki) can be 
done explicitly. Since we are considering the region where m 2 <C aim, an expansion in 
m, up to linear terms in m, can be used. The terms of order zero in m yield alternating 
signs when passing from one interval to the other and globally do not give significant 
contributions to the asymptotic domain. The order one terms in m (multiplied with the 
external m factor) can be treated perturbatively with respect to the leading aim factors. 
One finds the following behavior: 



7I\ o O 
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M* ~ (nr(n - n ) + (2 - -)m 2 - - ln(— ) + M 2 n > n , (A.26) 
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where uq is the m-dependent lower bound of n to be used. We have also included in 
the logarithmic term a mass independent contribution coming from the combination 
— f^'(0) + g 2 (ki) in the integrand. The coefficient of the m 2 ln(n/n ) term is only ap- 
proximate. 

We next consider the second region, where ann m 2 . This region is of relevance when 
heavy quarks or nonrelativistic expansions are considered. The above approximation for 
tanh (fii(ki) can still be used, but now making in the result an expansion in 1/m. The 
leading term of the expansion cancels the contribution of the factor ki of the integrand and 
the remaining terms are suggestive of a separation of the type M n = m + m n . The next- 
to-leading term contributes with alternating signs and is globally negligible. Retaining in 
the remainder the first leading term and the next one, the final result is: 



/ an \2/3 16 / an \ , /3m 2 \ n , . . 

1+2 fe) + 9^M^)J- n£n »- (A - 27) 



where we have chosen for the upper bound of n the same Uq as in Eq. flA.26j) . to ensure 
continuity between the two regions. The coefficients of the last two terms are approximate. 
The expression of the second term in the right-hand side of the above equation is typical 
of a nonrelativistic semiclassical energy of a particle of mass m in the presence of a 
linear confining potential (with a different coefficient), although we do not have in the 
present problem explicit external forces. The latter result at least justifies the validity 
of nonrelativistic expansions in domains where the number n respects the upper bound 
inequality imposed by the value of m. The mass dependence of the nonrelativistic term can 
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also be obtained by a direct analysis of Eqs. l lA.ip -l lA~~4"|) with the use of the corresponding 
expressions of the residues b n . 

The domain corresponding to the interval [0, Mi] does not belong to the asymptotic 
regions; nevertheless, it is useful to obtain for Mi an approximate explicit formula. Since 
that domain is concerned only by one mass, the corresponding approximation to be used 
for the functions /+ and /_ will be different. Here, one should keep in /_ its pole term 
in Mi and in / + its pole term in —Mi. No further approximation is needed, since 
in this case the function tanh<£>j(/cj) correctly vanishes at fc, = 0. Explicitly, one has 
ta,nh (fii(ki) = —ki/M\ and its integral can be evaluated. Keeping in Eq. ( 15.161) the factor 
[hi + mtanh(/)j(/cj)) as the integrand, one obtains for Mi the formula 

1 



Ml = - [m 2 + an + my/m 2 + 2an J, (A.28) 



which only provides the order of magnitude of Mi, but could be used in Eq. (1A.26j) when 
no = 1 (small or null quark mass). 

The asymptotic expansions ( 1A.26I) and ( 1 A. 2 71) were verified in numerical calculations. 



When m = 20Ja/7T, the value no = 30 seems to be optimal. For the other masses that are 



considered, m = 0., 0.1, 1., 5., in unit of ya/ir, n can be set equal to 1 and expansion 
1IA.26P can be used with formuula (IA.28|) for M 1 . 

A. 3 Numerical resolution 

The reconstitution of the functions F\ and F [Eqs. 1 13 .2p and 1I3.3P ] and of the functions 
sm(p and cosy? [Eqs. 1 15. 24|) and 1 15. 25j) ] necessitates the knowledge of at least a few 
hundred mass parameters M n . Several difficulties, however, are met during the numerical 
resolution of Eqs. llA^ - l[Pl) . 

Equations (1A. 1|) - (1A.4[) represent two sets of compatible equations concerning the same 
unknown quantities (essentially the M n s). Their compatibility was established when the 
whole contributions of the series were taken into account. Their truncation, for computa- 
tional purposes, introduces a small amount of incompatibility. Furthermore, each of the 
series has a different rate of convergence under iterative procedures. It is therefore not 
possible to solve separately and simultaneously these equations without meeting instabil- 
ity problems. The best way is to solve one of the compatible equations, or a combination 
of the two, and then check the numerical validity of the other. 

A second difficulty arises from the fact that the distance between two successive poles 



in the series decreases with increasing q, becoming of the order of yair/q. This means 
that the corresponding function varies very rapidly between two successive poles and small 
changes in M n induce big ones in the equations. 
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To analyze more quantitatively the accuracy of the calculation, we designate by e = 
a generic equation of the type of Eqs. (lA.ip and ( 1A.3|) . by e' = its derivative equation, 
of the type of Eqs. ( 1A.2|) and ( 1A.4[) . and by x one of the masses M n . Using Newton's 
method of iteration [25J, equation e' = leads to the following iterative relation: x^ 1 ^ = 



where e" is the second-order derivative of the function e with respect to x (M n ) 
and i is the order of iteration, e" is negative and for large n its modulus reaches values 
of the order of 10 4 — 10 6 (in unit of Ja/n). In these domains, the rate of convergence of 



the iteration is very slow and violations of the equation e! = by factors of unity or 10 
actually correspond to very small departures of M n from its exact value. The analysis of 
the equation e = is more complicated, since here the first derivative of e is itself null or 
small and one needs to continue the expansion up to the second-order derivative. 

A third difficulty arises from the weak convergence of the series of Eqs. (1A.1[) - (1A.4[) . 
It requires a detailed treatment of the asymptotic tails. 



For the numerical resolution, we have used the combination e—e'^a/nof the equations 
e = and e! = and checked after each iteration the separate validity of each of them. We 
have considered 900 active points M n and similarly for the m n s determined simultaneously 
through Eqs. (13.141) . As input values, the asymptotic expansions ()A.26j) and ()A.27j) are 
used. For the asymptotic regions of the various series, 700 background or passive points, 
from 900 to 1600, are considered using the asymptotic formulas (1A.26|) and (1A.27|) with 
a global scaling factor adjusted after each iteration by continuity at the frontier point 
n = 900. Beyond that region, the contribution of the series is approximated by an 
integral. High precision on the results is reached after a number of iterations that varies 
between 10 4 and 5 x 10 4 . The accuracy of the results concerning the first 400 points 
is better than 10 -4 and slightly decreases for the remaining points. The uncertainties 
are mainly due to the contributions of the asymptotic tails. Finally, various consistency 
relations have been checked: equations (13.131) . ( 1A.10I) . ( 1A.9I) (for negative values of z), 
and the constraint sin 2 (p(k) + cos 2 <p(k) = 1. 
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